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\ We have calculated the mass accumulation efficiency during helium shell flashes to 

examine whether or not a carbon-oxygen white dwarf (C+0 WD) grows up to the 
Chandrasekhar mass limit to ignite a Type la supernova (SN la) explosion. It has been 

■ frequently argued that luminous super-soft X-ray sources (SSSs) and symbiotic stars 
are progenitors of SNe la. In such systems, a C+0 WD accretes hydrogen-rich matter 

■ from a companion and burns hydrogen steadily on its surface. The WD develops a 
helium layer underneath the hydrogen-rich envelope and undergoes periodic helium 
shell flashes. Using OPAL opacity, we have reanalyzed a full cycle of helium shell 
flashes on a I.SMq C+0 WD and confirmed that the helium envelope of the WD 
expands to blow a strong wind. A part of the accumulated matter is lost by the wind. 

The mass accumulation efficiency in helium shell flashes is estimated as 
me = -0.175 (logM + 5.35) + 1.05, for -7.3 < log M < -5.9, and r? H e = 1 for 
—5.9 < logM < — 5, where the mass accretion rate M is in units of M & yr -1 . 
In relatively high mass accretion rates as expected in recent SN la progenitor 
models, the mass accumulation efficiency is large enough for C+0 WDs to grow 
to the Chandrasekhar mass, i.e., r/He = 0.9 for logM = —6.3, and r/He = 0.57 for 
logM = —7.0. The wind velocity (~ 1000 km s _1 ) is much faster than the orbital 
velocity of the binary (< 300 km s _1 ) and therefore, the wind cannot be accelerated 
further by the companion's motion. We suggest observational counterparts of helium 
shell flashes in relation to long term variations in super-soft X-ray fluxes of SSSs and 
symbiotic stars. 
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Subject headings: binaries: close — novae, cataclysmic variables — supernova: general 
— stars: mass-loss — white dwarfs 



1. INTRODUCTION 

Type la supernovae (SNe la) are widely believed to be a thermonuclear explosion event of 
a white dwarf (WD), although the immediate progenitors have not been identified yet. Among 
many evolutionary paths toward SNe la proposed so far, arguments are now focused on the 
Chandrasekhar (Ch) mass model vs. the sub-Chandrasekhar (sub-Ch) mass model, and the single 
degenerate (SD) model vs. the double degenerate (DD) model (e.g., [Branch et al. 1995 for a 



review of SN la progenitors). Of these models, the SD & Ch model is the most promising one 



both from the observational and theoretical points of view (e.g., Nomoto et al. 1994 for a review; 



Kobayashi et al. 1995 



One example of the SD & Ch models is luminous super-soft X-ray sources (SSSs) in which 
the WD accretes hydrogen-rich material at relatively high rates (M > 1 x 10 yr _1 ) (van 

den Heuvel et al, 1992). The WD accretes almost pure helium (Y = 0.98) as a results of stable 
hydrogen shell-burning and experience helium shell flashes. In a strong flash, the envelope expands 
greatly (e.g., Iben 1982^ ) and may suffer wind mass-loss or a Roche lobe overflow (e.g., Kato, et 



al. 1989] ). If a large part of the helium layer will be lost by the wind, the WD hardly grows in 



mass to the Chandrasekhar mass flCassisi et al 199~8 ). Thus, it is essentially important to estimate 



the efficiency of the net mass accumulation onto the C+O core. Kato, Saio, & Hachisu (1989, 
hereafter KSH89) has first estimated the mass accumulation efficiency in helium shell flashes. 
They used the old opacity, however, and it should be recalculated with the new opacity. 

A new evolutionary path of the SD & Ch scenario has been proposed by Hachisu, Kato, & 
Nomoto (1996; hereafter HKN96), in which a WD accretes hydrogen-rich matter from a lobe-filling 
red-giant (WD+RG system). In this model, optically thick winds blow from WDs, which stabilize 
the mass transfer even if the donor has a deep convective envelope. The WD steadily accretes 
hydrogen-rich matter from the companion. Its mass eventually reaches the Chandrasekhar mass 
limit and the WD explodes as an SN la. Hachisu, Kato, & Nomoto (1999a, hereafter HKN99) 
show this WD+RG system is one of the main channels to SNe la. Li & van den Heuvel (1997) 
extended HKN96's model to another type of binary systems consisting of a mass-accreting WD and 
a lobe-filling, more massive, somewhat evolved main-sequence star (WD+MS system). Hachisu &: 
Kato (1999) further proposed an evolutional path to binaries consisting of a WD and a somewhat 
evolved MS star characterized by helium-rich accretion as observed in U Sco. This is another main 
channel to SNe la as shown by Hachisu, Kato, Nomoto, &; Umeda (1999b, hereafter HKNU99). 

All of these models will experience helium shell flashes followed by optically thick winds 
during its way to an SN la. If most of the helium matter is blown off in the wind, the WD 
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cannot grow to the Chandrasekhar mass (Mch = 1.38Mq). Therefore, once again it is essentially 
important to evaluate the efficiency of mass accumulation during the wind phase of helium shell 
flashes in order to reach the definite conclusion of SN la progenitors. 

Cassisi, Iben & Tornambe (1998, hereafter CIT98) criticized HKN96's model and concluded 
that the white dwarf mass hardly increases to Mch because of heavy mass loss owing to a Roche 
lobe overflow during helium shell flashes. However, once the optically thick wind occurs during 
helium shell flashes, the wind velocity is as fast as ~ 1000 km s _1 , much faster than the orbital 
velocity (< 300 km s _1 ) so that the companion motion hardly affects the mass ejection. The 
wind matter quickly goes away from the system without any interaction to the orbital motion. 
Therefore their argument has no physical meaning when a strong wind occurs. 

In this Letter, we have followed a full cycle of helium shell flashes including the mass loss 
owing to optically thick winds in order to obtain the mass accumulation efficiency of accreting 
WDs. We have used the same techniques as in KSH89 but adopted OPAL opacity that has a 
strong peak around logT (K) ~ 5.2, which drives a strong wind as shown by Kato & Hachisu 
(1994). In §2, we describe our methods, assumptions, and main results of helium shell flashes on a 
1.3Mq WD. Discussions follow in §3. 



2. MASS LOSS DURING HELIUM SHELL FLASHES 

Kato, Saio, &: Hachisu (1989, KSH89) adopted two different approaches to follow a full cycle 
of helium shell flashes: time-dependent hydrostatic calculation and mimicking time-dependent 
behavior by a sequence of steady-state wind solutions. They combined these two methods, because 
the time-dependent Henyey code failed when the envelope expands greatly while the steady-state 
approach works well. We have used the same basic approaches as in KSH89 to follow a full cycle 
of helium shell flashes on a 1.3Mq WD. The main difference is in the wind solutions, which we 



replaced by the new solutions recalculated with the updated OPAL opacity ([Iglesias fc Rogers 



1996|) . 



The ignition models are unchanged because OPAL opacity gives essentially the same values 
as the old opacity does at higher temperatures of log T (K) > 6. We have assumed the same 
white dwarf radius of 1.3M0, i.e., logR (cm) = 8.513, and the envelope chemical composition 
of Y = 0.48, C + O = 0.5 and Z = 0.02 for helium, carbon and oxygen, and heavy elements 
(including C+O in solar ratio) by weight, respectively. Such a small helium content is expected 
as follows: in the early stage of a shell flash, the convection spreads over the entire envelope, then 
descends and disappears before the photosphere reaches the maximum expansion. The ashes of 
helium burning, i.e., carbon and oxygen, is mixed into the entire envelope to reduce the helium 
content by about a half on average from the initial value of Y = 0.98. Thus, we have assumed the 
uniform composition with Y = 0.48 in the quasi-evolutionary sequence. The convective energy 
transport is included in both the static and the wind solutions by the mixing length theory with 



- 4 - 



the mixing-length parameter of 1.5. The convection is, however, inefficient in the wind solutions, 
because a wind is too fast (supersonic) compared with the convective motion (subsonic). (See 



Kato & Hachisu 1994 for details.) 



Figure 1 depicts evolutionary tracks of one cycle of helium shell flashes. After the onset of 
a shell flash, the helium layer expands to reach the maximum expansion of the photosphere, the 
radius of which depends on the ignition mass AM; g of the helium layer. We call this rising phase. 
In the rising phase, the optically thick wind occurs when the photospheric temperature decreases 
to logTph ~ 5.45. After the maximum expansion, the photospheric radius decreases with the 
envelope mass and the star moves blueward in the HR diagram. After the wind stops, the envelope 
mass further decreases owing to nuclear burning and the star moves downward to come back to 
the original position. We call this decay phase. 

For the rising phase, evolutionary tracks are approximated by a sequence of static solutions 
with a constant envelope mass until the optically thick wind occurs (see KSH89 for details). Three 
tracks have been calculated with (a) AM ig = 4.8 x 1CT 5 M , (b) AM ig = 8.7 x 1CT 5 M , and (c) 
AM ig = 3.5 x 1O~ 4 M . These ignition masses are realized for the helium accretion rates of (a) 
M acc = 1.7 x 1(T 6 M yr" 1 , (b) M acc = 6.7 x 1O~ 7 M yr" 1 , and (c) M acc = 1 x 1(T 7 M vr" 1 , 
respectively. The helium accretion rates of (a) and (b) correspond to the largest rates of hydrogen 
steady shell-burning, M acc ~ LEdd/Xen, in the envelope with X = 0.35 and X = 0.7, respectively, 
where -LEdd on 1/(1 + X) is the Eddington luminosity for electron scattering and eh is the nuclear 
energy release of hydrogen-burning per unit mass. Thus M acc is proportional to 1/^(1 + X) for 
the hydrogen content X. Such a low hydrogen content of X = 0.35 (or Y = 0.63) is expected in 
the helium-rich matter accretion in recurrent nova U Sco ( [Hachisu fc Kato 1999| ). 

The evolutionary track for (a) AM- lg = 4.8 x 1O~ 5 M is very close to the curve for the decay 
phase (thick curve) and omitted in Figure 1. No optically thick winds occur in this case because 
the photospheric temperature does not decrease to logT p h ~ 5.45. At the maximum expansion, 
the photospheric temperature reaches logT p h = 5.83 (indicated by arrow a), and the radius 
-Rph = O.O22i? . 

In the case of (b) AM- 1S = 8.7 x 1O _5 M , the wind mass-loss begins at logTph = 5.45. The 
maximum expansion is reached at R p h < O.7Oi? and logT p h > 5.05 (indicated by arrow b), where 
the equality stands when we neglect the wind mass loss in the rising phase. The star moves 
blueward along the thick curve in the decay phase. The envelope mass is decreasing owing to 
wind and nuclear burning. The optically thick wind ceases at logT p h = 5.45. In the third case of 
(c) AM ig = 3.5 x 1O _4 M , the wind mass loss begins at logTph = 5.44, and the star reaches the 
maximum expansion at i? p h < 9.1i? and logTph > 4.45 (indicated by arrow c). 

These tracks for three different envelope masses merge into one track after the maximum 
expansion, because the envelope reaches a thermal equilibrium in the decay phase. The more 
massive the ignition mass is, the more redward the star moves in the HR diagram, that is, 
the star is traveling a longer track to a lower temperature. If the ignition mass is as large as 
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AM ig ~ 7.7 x 10~ 4 M Q , the star reaches logT ph > 4.07. This envelope mass, AM ig = 7.7 x 1O _4 M 
corresponds to the ignition mass at M acc = 4.5 x 10 _8 M Q yr -1 . 

Figure 2 shows the photospheric temperature T p h, the photospheric radius R p h, and the 
photospheric wind velocity V p h in the decay phase of helium shell flashes. The arrows indicate 
the same three ignition masses as in Figure 1. At the maximum expansion, the star reaches these 
position and then moves blueward as the envelope mass decreases owing to wind mass-loss and 
nuclear burning. The wind massdoss stops at AM = 5.0 x 10~ 5 M Q , marked by a filled circle. The 
photospheric velocity is as large as ~ 1000 km s _1 , which is much faster than the previous results 
(~ 300 km s _1 in KSH89), because of a strong acceleration driven by the large peak of OPAL 
opacity. 

The newly obtained wind velocities are much faster than the orbital velocities of SN la 
progenitors of the WD+MS or WD+RG systems (< 300 km s _1 ). Therefore, we can neglect 
effects of the binary motion on the mass ejection during the common envelope phase, because the 
envelope matter quickly goes away with little interaction with the orbital motion. This implies 
that neither a Roche lobe overflow nor a common envelope action plays an important role on mass 
ejection. 

Figure 3 shows the total mass decreasing rate of the envelope log \d(AM)/dt\, which is the 
sum of the wind mass loss rate M wim } and the mass decreasing rate owing to nuclear burning 
M nuc against the envelope mass log AM for the decay phase of helium shell flashes. The processed 
matter accumulates at a rate higher than the mass ejection rate for the region of log AM < —3.74, 
and then a large part of the envelope matter accumulates on the white dwarf. These results are 
very different from the case of nova outbursts, in which most envelope mass is blown off. It is 
because the difference is in the nuclear energy release per unit mass between hydrogen-burning and 
helium-burning: hydrogen burning produces energy ten times larger than that of helium burning, 
i.e., only a tenth of the hydrogen-rich envelope mass is enough to blow the entire envelope off. 

Figure 4 shows the accumulation efficiency defined by the ratio of the processed matter 
remaining after one cycle of helium shell flash to the ignition mass. The amount of matter lost by 
the wind in the decay phase is calculated from the rate of wind mass-loss and nuclear burning in 
Figure 3. This accumulation ratio is well fitted by the solid curve expressed by 

f 1, for -5.9<logM< -5 

VHe = < / • \ 2 . ~ f 1) 

[ -0.175 (log M + 5.35J +1.05, for - 7.3 < log M < -5.9 
where M is in units of M Q yr -1 . 

3. DISCUSSIONS 

Cassisi et al (1998, CIT98) have followed the evolution of hydrogen-accreting white dwarfs. 
These WDs are undergoing many cycles of hydrogen shell flashes until a strong helium shell flash 
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develops to expand the envelope to a red-giant dimension. Based on these calculations with the 
old opacity, they concluded that a hydrogen-accreting white dwarf will not become an SN la, 
because almost all of the helium envelope is lost from the system by a Roche lobe overflow or a 
common envelope ejection. CIT98 further denied the possibility of HKN96's scenario. 

In what follows, we criticize CIT98's argument and support the growth of C+O WDs as 
expected in HKN96's and Li <fe van den Heuvel's (1997) models. CIT98 calculated only much less 
massive white dwarfs of O.516M and O.8M to criticize ~ 1.3M Q model of HKN96. The trends 
of the envelope expansion and the mass accumulation depend on the WD mass because of the 
large difference in the surface gravity; nuclear energy release in helium burning is comparable 
to the gravitational potential energy for massive WDs, i.e., the ratio of nuclear /gravity is ~ 1.0 
for a 1.3M WD, whereas it is as large as 3.9 for a O.8M WD, and 8.9 for a O.5M WD. For 
massive WDs, therefore, a large part of helium is consumed to produce nuclear energy before the 
envelope expands against the gravity. It reduces the helium content of the envelope. This low 
helium content in the envelope increases the burning speed in mass coordinates as understood 
from M nuc ~ L^dd/Yeiie, where £He is the nuclear energy release per unit mass for helium burning. 
Thus, CIT98's criticism is inappropriate in the sense that their results for less massive WDs 
cannot be directly applied to very massive WDs. 

Once the optically thick wind occurs, a Roche lobe overflow or a common envelope ejection 
discussed in CIT98 becomes physically meaningless, because the wind velocity ~ 1000 km s _1 is 
much faster than the orbital velocity of the companion (< 300 km s _1 ) of typical SN la progenitors 
(WD+MS or WD+RG model). The wind goes away quickly without any effective acceleration by 
the companion. After the wind stops, the photospheric radius (< 0.1R & ) is much smaller than the 
Roche lobe (> IRq). Therefore, neither a Roche lobe overflow nor a common envelope ejection 
works in our binary system. 

We have ignored hydrogen-rich envelope above the helium layer. An accreting 1.3 M WD 
has a hydrogen-rich envelope of 2 x 10 -7 — 3 x 1O _6 M above the helium layer, depending on 
the mass transfer rate, chemical composition, and the other details of binary models. This is 
much smaller (< 0.01) than the mass of the helium layer at ignition and we have neglected its 
contribution to the mass accumulation efficiency. The hydrogen-rich envelope will be quickly 
blown off at a very early stage of helium shell flashes. 

Helium shell flashes may be observed as transient super-soft X-ray sources. For high mass 
accretion rates as expected in SN la progenitors, shell flashes are so weak that the surface 
temperature never decrease to less than ~ 1 x 10 5 K. Therefore, it will be a bright super-soft 
X-ray phenomenon with a faint optical counterpart. For a 1.3M WD, the bright SSS stage lasts 
~ 1 — 3 yrs with recurrence period of a few tens to several tens of years. For more massive WDs 
just before explosions, it lasts several months with the recurrence period of several years. 
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Fig. 1. — Evolutionary tracks of one cycle of helium shell flash on a 1.3 M white dwarf. Two tracks 
of the rising phase are shown by thin solid curves: the ignition mass is (b) AMi g = 8.7 x 10~ 5 Mc 







(upper), and (c) 3.5 x 10 M© (lower). The maximum expansion of the photosphere is indicated 
by an arrow for three ignition masses of (a) AM ig = 4.8 x 10~ 5 M Q , (b) AM ig = 8.7 x 1O~ 5 M , and 
(c) AMig = 3.5 x 1O~ 4 M0. The decay phases for these three ignition masses are merged into one 
line (thick solid). On the decay phase, the optically thick wind ceases at the point marked by the 
cross. The decay phase is extended to logT = 4.07 (short-dashed) for a strong helium shell flash 
with a large initial envelope mass AM; g = 7.7 x 10" 4 M Q . 
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Fig. 2. — Photospheric radius i? p h (cm), temperature T p h (K), and wind velocity Vp^ (cm s _1 ), 
are plotted against the envelope mass AM for the decay phase of helium shell flashes. Time runs 
from right to left. The wind mass-loss stops at AM = 5.0 x 1CT 5 Mq denoted by the dots. The 
three ignition masses shown in Fig. 1 are indicated by the arrows. 
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Fig. 3. — Wind mass-loss rate M w i n d (dashed), the mass decreasing rate owing to helium nuclear 
burning M nuc (dot-dashed), and the total envelope mass decreasing rate M w i n d + M nuc (solid) are 
plotted against the envelope mass AM. The ignition mass of three cases are shown by the arrows. 
Time runs from right to left. 
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Fig. 4. — Mass accumulation efficiency ryne, the ratio of the mass accreted after one cycle of helium 
shell flash to the ignition mass, is plotted against the helium accretion rate (filled circles). The solid 
curve is an empirical formula given by equation (1). The three mass accretion rates are indicated 
by the arrows, which are corresponding to the three ignition masses shown in Fig. 1. 



